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Abstract. This work is a further development of an approach to the description 
of relaxation processes in complex systems on the basis of the p-adic analysis. We 
show that three types of relaxation fitted into the Kohlrausch- Williams-Watts law, 
the power decay law, or the logarithmic decay law, are similar random processes. 
Inherently, these processes are ultrametric and are described by the p-adic master 
equation. The physical meaning of this equation is explained in terms of a random 
walk constrained by a hierarchical energy landscape. We also discuss relations 
between the relaxation kinetics and the energy landscapes. 
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1. Introduction. 

One of the distinctive features of complex systems (such as glasses, macromolecules, 
and proteins) is an anomalously slow non-exponential relaxation observed in a very 
wide range of time scales — from that of short-time molecular vibrations to days 
and months. Relaxation kinetics in complex systems is often fitted into three 
empirical laws, namely, the stretched exponential curve (the Kohlrausch- Williams- 
Watts law), ~ exp [— {t /t)"], < a < 1, the power decay law, ~ {t /r)~", < a, and 
the logarithmic decay law, a[ln {t /t)]^^ , 1 < a 0, ||, ^. These types of relaxation 
are characteristic of complex systems, just as exponential relaxation is typical for gases 
and liquids. 

Although formulation of a microscopic theory of anomalous relaxation in complex 
systems remains an issue of the day, there exist important concepts regarding 
the general relaxation behavior of complex systems. The starting point of these 
concepts is that the states of a micro-component of a system are strongly coupled to 
configurational rearrangements of the local environment, and it is important that these 
rearrangements cover a wide range of scales - from microscopic to mesoscopic. In order 
to describe configurational dynamics of this kind, the concept of energy landscapes 
is apphed According to this concept, configurational rearrangements are 

described by a random walk in the space of configurational states, and this process, 
in its turn, is determined by the energy landscape of the system. 

There is common understanding of the fact that a comprehensive description 
of energy landscapes of complex systems is hardly possible for several reasons. 
First of all, due to an uncertainty while selecting local potentials of interaction 
between the system's components. Another reason is the large number of metastable 
configurations, which occur due to multiple quenched constraints and the random 
nature of the local arrangement. The number of metastable configurations grows 
exponentially with respect to the size of the system. Thus, energy landscapes of 
mesoscopic systems (these systems are of the main interest) form multidimensional 
surfaces with an extremely large number of local minima. It is not quite clear 
how to describe those surfaces analytically and, accordingly, how to describe the 
configurational dynamics constrained by such landscapes. In this connection, 
computer simulation of model systems remains the main field of activity. These studies 
show that energy landscapes of complex systems may vary significantly. However, 
the characteristic types of relaxation obviously suggest that such landscapes have 
some common features. These features might turn out to be as fundamental for the 
dynamics of complex systems as the space symmetry is for crystal structures. 

One of the nontrivial assumptions, which initially appeared in the theory of 
spin glasses - [Q, was that the energy landscapes of complex systems have a 
hierarchical structure. In recent years, this conjecture was confirmed by the results 
of numerous computer simulation studies (see, for example, ^). It turned out 
that energy landscapes of complex systems represent hierarchically nested basins of 
local minima, viz., large basins consist of smaller basins of minima, the latter consist 
of even smaller ones, etc. The basins of minima are separated from one another 
by a hierarchy of activation barriers, i.e., the larger basins are separated by higher 
activation barriers, while the smaller basins within the larger ones are separated by 
smaller activation barriers, respectively. 

The existence of the basins of local minima allows one to describe configurational 
dynamics in terms of the basin-to-basin kinetics M . In this connection, it is assumed 
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that transitions between the basins belonging to the same level of the basin's hierarchy 
correspond to configurational rearrangements of a certain scale. Within such basins, 
a quasi-equilibrium distribution is established during relatively short (for a given 
scale) period of time, while transitions between the basins determine long-time 
configurational rearrangements. The hierarchy of nested basins corresponds to the 
hierarchy of scales of configurational rearrangements. 

As indicated in the early works on the spin-glass theory, the hierarchy of nested 



basins of local minima is similar to the hierarchy of nested ultrametric spheres |12 

. This idea can be found in some publications about anomalous relaxation in glass- 
like systems However, these investigations give no answer to the question 
whether ultrametricity can be directly used for the description of the dynamics of 
complex systems. 

Recent studies ||2^, ^ show that there is an approach which may lead to a 
noticeable progress in this direction. This approach is based on the p-adic analysis 
natural for ultrametric spaces. An introduction to the p-adic analysis can be found, 



for instance, in |22 . Some applications of the p-adic analysis in theoretical physics 
and theoretical biology are described in ||2^-||2^. 

The p-adic master equation describing a random walk constrained by the 
hierarchical energy landscapes of a certain type (the so-called ultrametric diffusion) 
was constructed in |2^, In these works, we describe a general method of finding its 
solutions and examine some Cauchy problems pertaining to conformational dynamics 
and reactions in proteins, and we also find some solutions of these problems expressed 
by analytic formulas. It is an interesting fact that the p-adic pseudo-differential 
equation introduced earlier as a p-adic analogue of the diffusion equation (see p2|) 
coincides with the particular p-adic master equation obtained in |2l[ |. 

Regarding p-adic pseudo-differential equations and stochastic processes in non- 
Archimedean spaces, see p2[, p^-pgj. 

The present paper is aimed at a further development of the approach suggested 
in j2^, for the description of configurational dynamics of complex systems. We show 
that the Kohlrausch- Williams- Watts law, the power decay law, and the logarithmic 
decay law can be described by p-adic master equations of the same form, and therefore, 
the relaxation laws characteristic of complex systems actually reflect the same type of 
ultrametric random processes. 

The paper is structured as follows. In Section 2, we introduce a general form 
of the p-adic master equation describing random walk constrained by a hierarchical 
energy landscape, and we also define a p-adic model of relaxation. In Section 3, we 
consider three versions of such models and establish their correspondence with three 
types of relaxation in complex systems. 

2. Random walk in ultrametric space 

The general p-adic master equation describing a Markovian process of random walk 
in ultrametric space can be written as follows: 

df{x,t) 



at 



[w{x\y)f{y, t) - w(y\x)f(x, t)] dn{y), (1) 



where x is a p-adic number, dfi(x) is the Haar measure on the field of p-adic 
numbers Qp, and t is time, t S IR.+ . Equation (|l]) is the usual balance equation 
for transitions between system states. The only peculiarity is that the space 
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of states is described by the p-adic numbers, i.e. this space is uhrametric. 
The function f{x, t) : Qp x R+ is a probabihty density distribution: the 

integral Jg /(x, t)d^{x) is the probabihty of finding the system in a domain i? C Qp at 
the instant t. The function w {x\y) : Qp x Qp R_|_ is the probability of the transition 
from the state y to the state x per unit time. 

The transition from a state y to a state x can be perceived as overcoming the 
energy barrier separating these states. The structure of the relation between the 
transition probability w {x\y) and the parameters of the energy barrier depends, in 
general, on the mechanism of configurational rearrangements, transition pathways, 
their distribution, etc. Finding this expression on the basis of the microscopic theory 
is a nontrivial problem beyond the scope of this article. For the sake of definiteness, 
we will use the standard Arrhenius relation 

u;(x|y)^A(r)exp|-^^| (2) 

where U {x\y) is the height of the (effective) activation barrier for the transition 
from the state y to the state x, k is the Boltzmann constant, and T is temperature. 
Studies of thermal behavior of anomalous relaxation in complex systems show that 
this expression is quite acceptable for a fairly wide range of temperatures. 

Formula (|^) establishes a relation between the structure of the energy 
landscape U {x\y) and the transition function w{x\y). In particular, the condition 
w {x\y) = w {y\x) corresponds to a degenerate energy landscape. An important 
example of a degenerate landscape is a regular hierarchical landscape. In this case, on 
each hierarchical level, the basins of states split up into the same number of smaller 
basins; moreover, the activation barriers between the basins of the same hierarchical 
level have equal values. Graphically, this hierarchical landscape is represented as 
a regular tree with the same branching rules for all branching nodes. As shown 
in |po| , the activation barriers U {x\y) separating the states x and y depend, in 
this case, only on ultrametric distances between the states, \x — y\^ = p'^ , 7 £ Z, and 
therefore, regular energy landscapes are completely described by functions of the form 
U{\x — y|p). Thus, in the case of random walk constrained by a regular hierarchical 
landscape, the master equation (|l]) takes the form 

^^^= f Wi\x-y\^)[f{y,t)-f{x,t)]d^^{y) (3) 

where 

W{\x-y\^) = j^^e^p 

^ F-y|p 

The p-adic master equation admits analytic solutions for a fairly wide class of 
functions W^(|a; — y\p). The technique of finding its solutions is described, for instance, 
in(§,|2^. 

It may seem that regular hierarchical landscapes represent a very idealistic model 
of real landscapes and have little relation to the latter. Indeed, regular hierarchical 
landscapes reflect no details of the energy landscapes of complex systems other than 
their hierarchical structure. Oddly enough, this is sufficient for the description of the 
characteristic types of relaxation. 

On the basis of general physical considerations, one can identify three essentially 
different types of regular hierarchical landscapes: 
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"Logarithmic landscapes" are characterized by activation barriers having a 
slow (logarithmic) growth with respect to the number 7 of the hierarchical level, 
U (7) ~ In 7. This type of landscapes may be associated with "flexible" systems having 
a wide range of available scales for configurational rearrangements. 

Contrary to logarithmic landscapes, "exponential landscapes" are characterized 
by a rapid (exponential) growth of activation barriers, J7(7) ~ e'* . Landscapes of 
this type represent "stiff" systems with a small set of distinct scales available for 
configurational rearrangements. 

Finally, there are "linear landscapes" with linear growth of activation barriers, 
U (7) ^ 7. These landscapes may be associated with systems of an "intermediate" 
type. 

In all these cases, we consider the Cauchy problem for the respective master 
equation (|^) with the initial conditions: 



A relaxation process will be understood as the evolution of population in the domain 
of the initial distribution \x\ < 1: 



In each of the above cases, the analytical estimation of relaxation function S{t) is being 
found. It would enable the relation of S{t) to three empirical laws for characteristic 
relaxation. 

3. Characteristic relaxation types 




(4) 




(5) 



The solution of the Cauchy problem for equation 
the form 



(||) with the initial condition (|^) has 



00 




(6) 



n=0 



where 




Indeed, applying the p-adic Fourier transformation to (||) and taking into account the 
initial condition (Q), we get 




'p 



where 




(7) 



Calculating the integral and using the substitution = p ^, we obtain 



00 



VF(p-") = (1 -p"i) p"W^(p'"+"+^) +P"^T^(p"+^) 



m=0 
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and therefore, (g) is the desired solution. 

In order to find S{t), the solution (^) should be substituted into (||). Calculating 
the corresponding integral, we find that 

oo 

5W = (p-l)5]p-"e-^"*. (8) 

n=l 

Here, for the sake of convenience, we have introduced the positive function 

oo 

A„=p" {l-p-')Y,p^"W{p"^+")+p-'W{pn . (9) 

m—O 

Next, we obtain estimates of S{t) for the above three regular hierarchical 
landscapes. The transition probability W is defined here by (^. The frequency 
factor is taken equal to the unity, and the temperature dependence is characterized 
by the parameter a = Tq /T . 

3.1. Logarithmic landscape. 

In view of (2), the transition probability function for the logarithmic landscape can 
be represented as 

1 1 

P \x-y\p In (1 + 

In this case, 

oo 

A„ = (1 - p-i) ln~"(l + p-^ ln-"(l + p") . (10) 

m=0 

Then, the following estimates hold: 

„, ln~"p 1 a\iC" p 1 



(a - 1) (n+ l)°-i (a-l)?i°-i 

Using these inequalities, we can estimate the relaxation function S{t) as follows: 

^\ 2a / t \^\ 
I a— l\lnp/ I 

<5(„<p^(exp{--l^(J-)*}+cxp{-(ji- "^^^ 

Thus, the long-time behavior of relaxation kinetics constrained by logarithmic 
landscapes is fitted into the Kohlrausch- Williams- Watts law. 

It is interesting to note, that as a — *■ 1, the upper bound for S{t) approaches an 
exponent. For a < 1, the series ( |l0| ) is divergent. This fact may be interpreted as a 
transition to exponential relaxation when a critical temperature Tq is attained. This 
kind of critical behavior of ultrametric diffusion constrained by logarithmic landscape 
was previously predicted in p6| on the basis of some qualitative considerations. 



p-adic description of characteristic relaxation in complex systems 7 

3.2. Linear landscape. 

In the case of a linear landscape, the transition probability function has the form: 

Estimates for the probability density distribution f{x,y) in the case of a linear 
landscape were obtained in ||2l|. Calculations for S(t) similar to those of yield the 
estimate 

ir (^i + 1^ {-Tp{^a)t)-i < S{t) < r (^i + 1^ {-Tp{-a)t)-i 

where r(a) is the gamma-function and ^p{a) is the p-adic gamma-function (see p^). 

Thus, the long-time behavior of relaxation kinetics constrained by linear landscapes 
is fitted into the power decay law. 

3.3. Exponential landscape. 

In this case, the transition probability function is 

W{\x~y\ )^y^e-^\'-'-y\^ , a>0, 

and from we have 



The following estimates hold: 



exp [-ap"] < A„ < _ ^ exp [-ap"] . 
Hence, for relaxation function S{t) we get 

. , "// ^ < S(t) <{v^ +p + e)-^ . 

2e Ini + ln(p"/(p" - 1)) w f ) 

The last estimate shows that the long-time behavior of relaxation kinetics constrained 
by exponential landscapes is fitted into the logarithmic decay law. 



It has been shown that the Kohlrausch- Williams- Watts law, the power decay and the 
logarithmic decay laws reflect the same type of random processes. Such processes 
are ultrametric and can be described by the p-adic master equation formulated 
in It seems that ultrametricity is a most general and fundamental feature 

of the dynamics of complex systems: The configurational space of such systems may 
be totally disconnected due to numerous quenched constraints, and therefore, it is 
impossible to accomplish a relatively large configurational rearrangement step-by-step, 
by way of small configurational changes. 
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In conclusion, we consider it necessary to make a few remarks. 

A direct comparison of the function S{t) with experimentally observed relaxation 
kinetics requires a certain caution. Indeed, the function S{t) is the population of 
a region in the configurational space, whereas in most experiments, the focus of 
observation is some probe coupled to configurational rearrangements, rather than 
configurational rearrangements themselves. A specific form of the relation between the 
probe state and a local configuration may happen to be important for the description 
of experimental data. 

The above three characteristic relaxation laws do not exhaust the diversity of 
relaxation kinetics in complex systems. There are numerous observations of "mixed" 
types of kinetics, with one type of relaxation kinetics observed in some time- window 
and another type observed outside. For instance, it may happen that at low 
temperatures the observed kinetics is well fitted into the logarithmic decay law, and at 
higher temperatures it is described by the power decay law. In other cases, the power 
decay may transform into the Kohlrausch- Williams- Watts law. It is not difficult to 
imagine energy landscapes whose activation barriers near the ground state are close 
to the exponential landscapes, while higher activation barriers correspond to a linear 
landscape and then transform into a logarithmic landscape. The above approach can 
be easily used for the description of ultrametric diffusion constrained by such " mixed" 
landscapes. Equation (|^) can be solved analytically for a wide class of functions 
W{\x — y\p , a) admitting the p-adic Fourier transformation. 

Finally, our considerations here are restricted to degenerate energy landscapes. 
These may be directly related only to some specific complex systems, for example, 
spin glasses. However, energy landscapes of other complex systems are obviously non- 
degenerate. Local minima in such landscapes are also clustered into a hierarchy of 
basins, and as the whole, they look like a global potential "hole" (or "funnel") with 
extremely rugged "walls" . The description of random walk on such landscapes by p- 
adic equations is of great interest, particularly, in connection with the applications of 
the concept of hierarchical energy landscapes to the protein folding |Q (see also | pl[ 
and references therein). Some results regarding this problem will be presented in our 
further publications. 
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